NACA TN4188 iggnj 


TECH UBRARY KAPB, NM 



NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 4188 


CHARTS RELATING THE COMPRESSIVE AND SHEAR BUCKLING 

STRESSES OP LONGITUDINALLY SUPPORTED PLATES 

TO THE EFFECTIVE DEPLECTIONAL STIFFNESS 

OF THE SUPPORTS 

By Aldle E. Johnson, Jr. 

Langley Aeronautical Laboratory 
Langley Field, Va. 



Washington 
February 1958 


L. - snARY 

AFL 2811 




N^ITIQMAL ADVISOEiy CCMCHraEE FCR AERQMAIJTICS 


TECHNICAL NOTE 4l88 


tech library fCAPB, NM 



□□LLfldB 


CHARTS RELATING THE COMPRESSIVE AND SHEAR BUCKLING 
STRESSES OF LONGirODINALLY SUF^TEH) FIA!1SS 
TO 'TTTB! EITECTIVE INFLECTIONAL STIFFNESS 
OF QHE SUPPORTS 
By Aldle E. Johnson, Jr. 


SUGARY 


A stability analysis Is made of long flat rectangular plates 
subjected to botb shear and CGn^>re8slve loads. The edges of the plates 
are taken to be Blxoply supported and the plates are supported along one 
or two Intezmedlate longitudinal U^es by lines of deflectl on al springs 
(elastic sv^orts). The results of the analysis are presented In 

the form of charts that are useful In the determination of the buckling 
load of plates stiffened by stringers or webs. 


INTRODUCTION 


A problem frequently encountered In the design of aircraft struc- 
tures Involves the determination of the shear and ccmpresslve stresses 
that a stiffened plate can suBted.n without buckling. Design Information 
Is available (refs. 1 to 5) on the stresses attelnable In a stiffened 
plate before the onset of buckling due to either a shear or a caznpresslve 
loading; but little Information Is available on the stresses attainable 
In plates subjected to combined shear and ccm^esslve loadings. Flat 
simply supported rectangular plates stiffened by one or more longitudi- 
nal stiffeners were analyzed In reference 1 for a shear loading and In 
reference 2 for a compressive loading. The present paper treats the 
CELse of oomblned shear and compressive loadings on long simply supported 
plates stiffened by one or two longitudinal stiffeners. 

In the analyses of references 1 and 2, stiffeners are Idealized as 
beam columns whose support to the plate can be e^ressed In terms of the 
area and moment of Inertia of the stiffeners and the wave length of the 
buckles. Practical construction, however, often Involves stiffening 
members whose support to the plate cannot be described by these parameters. 
Reference 3, therefore, presents a stability analysis of stiffened plates 
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In ccmpresslon In terms of a "generalized support stlffliesB parameter" 
vhlcli permits the Inclusion of cross-sectional and sheaz*lng distortions 
of the stiffening member In an evaluation of the support offered the 
plate. The stability criteria of the present Investigation are presented 
In terms of this generalized parameter. Ohe charts given are similar to 
those of reference 3 and extend portions of that work to include the 
euldltlon of shear loadings. 


SYMBOLS 


A, B 
A/bt 

^n 

b 


D 

E 

El 

bD 

K 


Nx 

n 


coefficients defining amplitude of support deflection 

stringer area ratio 

Fourier coefficients 

width of bay between support lines 

ratio of average stress In stringer to average stress In 
plate 

•r+3 

plate flexural stiffness per unit width, — 

12(l - ^2) 

Young's modulus of eleistlclty 
bending stiffness ratio 


coefficient 

nondlmenslonal conpresslve buckling -load coefficient, 


j£2d 


N_,b2 


nondlmenslonal shear buckling -load coefficient, 

coaopresslve load per unit width reg.ulred to cause buckling 
sheBo* load per unit width required to cause buckling 
Integer 
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Q 

t 

Ui, U2 

Vl, V2 
w 

ac, y 

P - S/-b 

Aj 

X 

tJL 

♦ 

+b5/«^ 


energy peuraxoeter 

thickness of plate 

strain energy of defomatlon 

work of external forces 
deflection normal to plane of plate 

coordinate axes In length and width directions of plate ^ 
respectively 

Leigranglan multipliers 

length of buckles 
Poisson’s ratio 

deflectlonal stiffness per unit length of support 
nondlmenslonal deflection restraint parameter 


STATEMENT OF PROBLEM 


The structural configurations analyzed herein are shown In figure 1. 
They represent long simply supported plates which are stiffened by one or 
two longitudinal stringers or full-depth webs and which are subjected to 
shear and compressive loeids. The stringers or webs (shown els springs) 
provide support to the plate which is simulated for the buckling mode 
considered herein by lines of springs wl-tti deflectlonal stiffness per 
unit length +. The mode of instability that Is of interest Is one in 
which the intermediate supports (lines of springs In fig. 1) deflect eua 
the plate buckles Into a series ^ skewed troughs and crests. 

The use of lines of deflectlonal springs to simulate the si;q?port 
provided by longitudinal stringers and full-depth webs Is discussed In 
reference 5 and utilized in refearences 5 to 5 "to predict the behavior of 
aircraft structures. Reference 3 also gives the effective spring stiff- 
ness of stringers and webs In terms of the buckle lengthy the stress 
level, and geometrical properties of those menibers. As an example, the 
effective spring stiffness of a longitudinal stringer with a sturdy cross 
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sectloa (and hence one whose deflectlonal stiffness Is defined hy ele- 
mentary beam theory) Is given by the following expression: 


»b3 ^ 1./EI . 




( 1 ) 


The nondlmeiLBlosal spring stiffness in equation (l) is given 

jc^D 

in terms of nondlmenaional parameters cananonly en^loyed in stability 
analyses of stiffened plates. One of these parameters^ is a function 

of the others appearing on the rig^t-hand side of the equation. (See 


axe available in references 3 


ref. 6.) Corresponding expressions for 

to 5 for cases in vhich the distortions of attachment flanges of supporting 
members play an Inrpgrtant role in determining their effective spring stiff- 
ness. The reader is referred to these references for detailed discussion 
of this and related problems in evaluation of support stiffnesses. For 
sinpliclty; equation (l) is used in discussing applications of the results 
presented herein. 


RESULTS 


Stability criteria for the structures depicted- in figure 1 are 
derived in the appendix by the Lagrangian multiplier method. crite- 

rion for the two-bay structure (fig. 1(a)) is given by equation (A52) 
and the criterion for the three-bay structure (fig. 1(b)) is given by 
equation (A28). T^ nondimensional parameters appearing in these criteria 
are the buckling coefficients kg and kg^ the buckle aspect ratio 


and the stiffness parameter 



The deflection functions used to derive the stability criteria are 
given by eqioations (Al) and (A29) and have symmetry about certain points 
(the crests and trou^^ of buckles) on the longitudinal center line of 
the plate. Solution of the stability criterion for the two-bay plate 
yields a single root corresponding to the first BynKnetric mode - a half- 
wave across the plate. Solution of the stability criterion for the 
three -bay plate yields two roots corresponding to the first two symmetric 
modes. However, only the results for the first symmetric mode are pre- 
sented because this mode occurs at a lower stress level than the second 
synanetric mode and is therefore more likely to occur unless additional 
restraints are introduced along the support lines. 
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A Numerlc 8 Ll calculations vex^ perfozmed on the National Bureau of 

>|fh^ 

Standards Eastern Automatic Ccxiiputer to determine the value of 

associated with assigned values of and kg. A typical set of 

data Is given In table I. 0!he data vere reduced so that^ for each con* 
figuration considered, design charts with one of the stress coefficients 
held constant could be constructed. ^Ehe charts give, for constant 
values of the support stiffness, the second stress coefficient sis a 
function of the buckle length at vhlch buckling occurs. The charts are 
presented In figures 2 to 


Two -Bay Plate 

are given In flg- 
of and p 

kg In terms of 
» 0 ) was taken 


slniply si^ported 

plates stiffened by a single longitudinal stringer can be obtained from 
the data of figure 2(a) (kg » 0) of the present report. The value of ^ 

given In reference 1 for a given stress level Is the maxlTnum value of 
^ obtained with the use of equation (l) and the combinations of 

and p of figure 2 (a) of the present report at the stress level under 
consideration . 


The results of calculations on the two -bey plate 
ures 2 and 5. Figure 2 gives values of kg In terms 
for k^ held constant, £uid figure 3 gives values of 
and p for kg held constant. Figure 5 (a) (kg 
from the data of reference 5 . 

The results given by figure 3 of reference 1 for 


Ihree-Bay Plate 

^[he results of ceilculatlons on the three -bay plate are given In 
figures 4 and 3« Figure 4 gives values of kg In te 2 mis of and 

for kg held constant, and figure 3 gives values of kg In terms of 


— and p for kg 
data of reference 3 . 


held constants Figure 3(^) vas taken from the 
As In reference 3 a cutoff to the curves Is 


P 
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Included In figure ^(a) because above this cutoff the dominant mode of 
Instability changes from one which Involves deflections of sig>port to 
one In which buckling takes place with no deflection of the siQ>x>orts 
(local buckling In each bay). The curves determined by the present 
analysis are shown as biToken lines above the cutoff In order to provide 
Informatl^ necessary for the solution of problens In the neighborhood 
of the cutoff. Ill a manner similar to that for the two-bay plate, the 
data of the present report (fig. 4(a) and eq.. (l)) can be used to obtain 
the curve for the three -bay plate given by figure 3 of reference 1. 

DISCUSSION 


The charts of figures 2 to 5 were designed to be applicable to 
plates long to accommodate several buckles but can be apx)Ued 

with reasonable &LCcuracy to problems Involving plates which buckle In 
as few as two or three buckles. Their use on shorter plates gives con- 
servative results. The degree of conservatism may be large or small 
depending upon the problem under consideration. In the case of rectan- 
gulEkr plates subjected to compressive loads, the plate Is te rm i n ated by 
ribs which contact the plate on a line that Is also a natural nodal line 
of the buckled plate, and the results obtained from the charts are exact. 
Plates stiffened by sturdy longitudinal stringers and chordwlse ribs 
usually develop buckles whose length Is large In comparison with the 
width of the plate as will be demonstrated later. The rib spacing for 
these structures may be such that the plate buckles Into a single half- 
wave between ribs, and use of the charts will give conservative results 
except In the case of a pure compressive loGidlng. Plates stiffened by 
stringers or webs which do not behave as sturdy stringers often develop 
buckles whose length is of the order of the plate width. The charts can 
be used for these structures with good results because ribs are usually 
several plate widths apart. The degree of conservatism may be large, 
however, for plates with «Tn«.n rib spaclngs which are subjected to shear 
stresses. For Instance, as an extreme example. If the shear buckling 
coefficient for a simply supported sq.uMe plate is obtained from the 

1-1 ^ \lfb^ 

charts presented herein |from fig. 2(a) with '' '' 


^rom fig. 

fig. 4(a) vlth 3 = 5 and = oV 
the values of k 3 read from figures 2(a) 


p B 2 and 


\|fb-^ 


0 or fr om 


a value of 5.5 is obtained etfter 


and 4(a) are multiplied by 2^ 


and 3^, respectively, to convert them to buckling coefficients In which 
the plate width b is the distance between the simply sxapported edges. 
The correct value for this case Is 9.3. The difference between these 
two numbers reflects the effect of restraint on the buckle shape near 
the ends of the plate. 
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p Aq application of the oharte presented herein Is given In figures 6 

and 7 where Intereictlon curves for two-bay and thiee-bay plates stiffened 

by sturdy longitudinal stringers of particular proportions are presented. 

% The curves were drawn with the use of flgilres 2 to 5 and eq.uation (l). 

The discontinuities In the Interaction curves result from a change In the 

mode of Instability analogous to the change In mode responsible for the 

cutoff In figure ^(a). [Qie iq>per section of a given curve ^curve for 

EI\ ' 

constant —1 corresponds to a mode of Instability In which the supports 

deflect with the plate. The lower section (from ref. 7) corresponds to a 
mode of Instability In which buckling takes place without deflection of 
the supports. Thus, from such curves, tlie mode of instability and the 
combination of shear and corngpresslve stresses that can be carried by an 
Infinitely long longltudlnal3y stiffened plate can be determined directly. 

Another application of the design charts presented herein Is given 
In figure 8 where the data of figure 2(a) ere presented In an alternate 
form which Illustrates the behavior of plates stiffened by sturdy longi- 
tudinal stringers. The results shown for the peirtlcular case of a two- 
bey plate loaded In shear are g.ualltatlvely similar to those for plates 
with two or more bays loaded In con^resslon or In coiiQpresslon and shear. 
When the bending stiffness of the longitudinal stringer Is zero (that Is, 

» the stringer Is Inactive or nonexistent), the plate buckles as a slnrply 

supported plate at k 3 ■ 1.J5 “od p = 2.4. The corresponding critical 

stress and buckle length check those found In other Investigations 
(ref. 8, for example). As the bending stiffness of the stringer Is 
Increased, the critical stress and wave length Increase. Finally, In 

this process a value of ^ can be found (** 1 , 750 ) so that the plate may 

bD 

buckle Into either of two modes with values of p that differ by a fac- 
tor of about 10 for the peorblcular case considered. For larger values 

of the plate buckles at the shorter wave length at values of kg 

which approximate those ccmputed for a plate simply si:pported at the 
stringer. In the conparable problem for a two-bay plate loaded In com- 
pression, the buckling coefficient for the deflection function used 
approximates that of a plate clamped at the stringer when the bending 
stiffness of the stringer Is large so that the plate bvickles at the 
shorter wave lengths. 

Further applications of charts like those presented herein are pre- 
sented In considerable detail In reference 5* 
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CONCLUDINa REMARKS 


Design charts have been presented which evaluate the camhlnation 
of shear and cciDpresBlve stresses required to buckle long flat simply 
supported rectangular plates stiffened by one or two Intennedlate lines 
of deflectlonal springs (elastic line supports). The charts are appli- 
cable to the detennlnatlon of the buckling loeid of plates stiffened by 
stringers or vebs of proportions encountered In aircraft construction. 


Langley Aeronautical Laboratory, 

National Advisozy Coomilttee for Aeronautics, 
Langley Field, Va., September JO, 1957 * 
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AEPEHSTDIX 

EERIVATICN OF STABILIT2r CRITERIA 


Stability criteria are derived by meems of tbe lAgranglan multiplier 
method for the buclclixig under combined ajdal coai^resBion and shear of 
iJoB simply supported plate configurations shown in figure 1. 


Three -Bay Plate 

The deflection function used to represent the buckled shape of the 
plate of figure 1(b) is given by the series 


w 



aln 2!Ey + cos 


sex 

7^ 


00 





b^cos^ 


(Al) 


with the origin of the coordinates as shqwn in figure 1. This deflection 
function can provide the desired buckle shape, but the req,uired geometric 
boundary conditions of zero deflection at the edges Eire not satisfied 
term by term. The required conditions, which can be written 


w(x,0) = w(x,3b) * 0 (A2) 

will be introduced by means of Lagrangian multipliers. 

The deflection of the first line of deflectional springs w^ can 
be represented by 

Wb = A sin 2^ + B COB ^ (AJ) 

A A 


From synnetry considerations, the deflection of the second line of 
deflectional springs Wg^ can be represented by 


V2b 


A sin 2^ - B cos 

A A 


Uk) 
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!Ihe condition that the deflection of the aprings is equal to the deflec- 
tion of the plate along the attaohznent lines 


w^(x) - v(x,b) 

(A5) 


(a6) 


will also he introduced by means of Lagrangian multipliers. 

The internal energy of the plate -support system is the sum of the 
strain energy of bending of the plate and the energy of the deflectional 
springs. The strain energy of a plate of stiffness D is 



(A7) 


where K is an undetermined coefficient that arises from the '^exact" 
differentiation of the series. The energy of the deflectional springs 
which have a spring constant if is 


The external work is in two parts. The work done by the shearing 
forces is 


Vn - -N. 


xy 


£ I t) 


dx dy 


4 ’^ 


t: 


an>n“ 


(A9) 
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The vork done by the uniform ccmpresBlve forces Is 


(AlO) 


Then 


+ Ua - Vi - Va) 


vhere 




• • • 


2: k 


• • • L 

V 9 y 9j 

5 


•N 

y 

GO. 

an^^ E V 


» • n»l,5>5j • • • 

Nxv-b^ 

H„b2 


(All) 


and 



Is a nondlTnens lonal deflectional spring stiffness pareineter. 


In order to satisfy the boundary condition of zero deflection 
(eq.. (A2)) and the condition that the deflection of the springs eq.uals 
the deflection of the plate (eq.s. (A 5 ) and (A 6 )), It Is necessary to 
Impose the following constraint relationships: 

bn = 0 (A12) 
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> an Bln ^ - A - 0 

n»l,5^ 5^ • • • ^ 


w 

5~ Tij, COB ^ - B - 0 

0^1,33,... ^ 


According to the Lagraoglan multiplier method, the esepression to 
he minimized 1 b 

-*)- 


j I: ^ 

V=l>3>5>*” / V=ij5i5> ••• / 


(A15) 


where ^ 2 ^ and axe the Lagranglan multipliers • 

0!he expresBions for Q must ndnlmized with respect to h^, 

K, A, B, A 3 _, ^ 2 j and A^« Minimization yields the following 

equations : 


(n “1, 5j • • • ®®) (A16) 




1 * 


2p3kgana 


Ag COB ^ + Aj = 0 (n * 1, 3, 5,, . . oo) (AI7) 
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2k 


r 

4-= 9 81 


K > 0 


0^,3,^,... 


3Q 4p^A I|fb5 


as 3 ,4 d “2 


§^“ a a^Bln^-A = 0 

••• 


o^l^3>5i ••• 


B » 0 


5%- 2_ ‘»-° 

^ n»l,5,5, ... 


(A18) 

(A19) 

(A20) 

(A21) 

(A22) 

(A23) 


EquatlooB (Al6), (AI 7 ); (A19); and (A20) may be solved for eIq, 
k, and respectively^ and the resulting egressions substituted 

Into equations (Al8), (A21), (A22), and (A23). This substitution results 
In the following set of simultaneous hom^eneous equations: 



A^Cn sin SSX, 

--L L 




cos 
- Cn^ 




00 

H 


\pA,3t^t • • 


AnBa 


Bn^ - Cn^' 


(A24) 
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(A25) 


(A26) 


(A27) 



Equations (A24) to (A 27 ), vben ccoiblned Into deternrlnantal fonij yield tbe foUovlrig 
stability crlterton: 



(MO) 


vhere 

► 

Two-Say Slate 

I 

followlzig deflection function Is used to represent the buckled shape of tbe plate of 
figure '1(a) : 



h- 
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T.R^r«ngiAn multlpUers are used to Introduce tlie cozidltlans of zero deflection at the edges 
of the plate azid the conditloQ. of con^tlhllity of deflection of plate and. springs a lon g the 
attachment line. !Qiese conditions are 


v(xjO) - v(x,2b) ■ 0 


(A30) 




(A5L) 


Internal eiiergieB and external vork are confuted In the same vay as for the three-bay 
plate and are. minimized vlth respect to the saiBe coeffidenta • Hie stability criterion thus 
detezndjaed is 
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Lf5f5j 
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(a) kQ = 0. 

Figure 2.- Stability curved for tvo-bey plate irlth constant coomreBBlve load. 
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(a) kg 0 . (Data are from ref. 3. ) 
n^ire 3.> Stability curves for tvo-bay plate vith constaat shear load. 
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I Figure 4.- Stability curves for three-bay ]^te with constant caqpresslve load. 
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Figure 4.*- Continued. 
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(a) kg a 0. (Data are from ref. J*) 

Figure 5.- StaMllty curves for three-bay plate with constant shear load- 
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Figure Concluded. 
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Figure 7-- ^^ypical Interaction curve for shear and con5>resslve buckling 
stresses of three -bay plate ’. A/bt =» 0.4. 
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